T 


i 


r 


&c>o £>/?'? ^ 

/ /j ' ' 6 


Semi-Annual Reports on 

Mesospheric Dynamics and Chemistry from SME Data 


by 


Darrell F. Strobel 
Principal Investigator 
Department of Earth and Planetary Sciences 
The Johns Hopkins University 
Baltimore, MD 21218 


July 1, 1986 to December 31, 1987 


NAG 5-796 


< MAS A-CB- 182465) HESOSPHEBIC ESMABICS AMD 
CEEMISTBI FECK SMI DATA Settiacnoal Report, 1 
jul. 1S6 1 - 31 Dec. 1S€7 (Jehus Hopkins 
Cniv.) 18 P CSCL 04A 


G3/46 


N88-* 17 139 


Dnclas 

0123405 



1 


The principal research efforts accomplished under NASA grant NAG 5-796 
consisted of two parts, which are described in the following two sections. 

1. Approximate Algorithm for 
Computing C0 2 I5p Band Cooling Rates 

A fast Curtis matrix calculation of cooling rates due to the 15^ band of 
C0 2 is modified to parameterize the detailed calculations by Dickinson (1984) 
of infrared cooling by C0 2 in the mesosphere and lower thermosphere. His 
calculations included separate NLTE treatment of the different 15^ bands 
likely to be important for cooling. 

Our goal was to compress the detailed properties of the different 15^ 
bands into a modified Curtis matrix, which represents one composite band with 
appropriate averaged radiative properties to allow for a simple and quick 
calculation of cooling rates given a temperature profile. 

We use the two-level approach of Houghton (1977) and write the ratio of 
the collisional deactivation rate of the upper level a 21 , to its spontaneous 
decay rate, A 21 , as 


<p= — = P/0.04mb 

and introduce the altitude dependent coefficient 

£* = ( 4nSn<t>y' 

with band strength S and concentration n at discrete heights, k. 

The Planck function B(T) , the black body emission per unit solid angle 
per unit area at temperature T, is related to the cooling rate Q, and the 
total NLTE source function J, by 


J k B k + E kQk 
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Q k 'y.C kjj i + Q SOLAR, k 
J 

where C kJ is the element of the Curtis matrix and Qsqlar, k *- s t * ie s °l a ^ heating 
rate (Houghton, 1977). At 15/i solar heating can be ignored. In matrix 
notation the latter equation may be written 

Q = CJ (1) 

but in the former equation there is no summation and it is incorrect to use 
matrix multiplication, unless we generalized non-diagonal terms. The source 
function becomes: 


J = 5 + b 


where 


(b) t -£*<?* 

We introduce the matrix Q i ^ with i the index for height level and j the index 
for latitude and define 


«<?„) 

J = B + b 

There is no easy way to invert matrix equation (1). Two methods of 
solution were tried: 

Method 1: The Q dependence of b is assumed to be slowly varying in T, 
with most temperature information contained in B. Thus b is held constant. 
This method is only valid when b- l db/dT<zB'dB/dT . 

Method 2: This is an iterative technique in which the temperature profile 
is used to generate B and the first estimate of Q is given by Q 0 = CB, J l = B 
+ b Q , etc. where B is the source function. For the next iteration the source 
function is J x = B + b Q (Q 0 ) and « CJ 1# This procedure is repeated until 
convergence is reached. 
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The following expressions were used 

5 = 97. 1 (s)^exp^^p I j E = (4nSntj>y' 

5 = 242 ./^ 

V 273 

J *rr 2 * KliE ’ Q)li 

with additional, adjustable constants introduced to achieve the best 
accuracy. We selected K t = 4.11 x 10 11 and K 2 * 81.12. These values were 
chosen because they kept the diagonal Curtis matrix elements near unity. 

Inverting the equations to obtain the Curtis matrix C yields an unstable 
matrix, such that the smallest change in the temperature profile produces 
large fluctuations in the cooling rates, clearly not desirable. 

Our initial trial guess for the Curtis matrix was the identity matrix 
minus the same constant for all terms, both on and off the diagonal. We 
assumed that half the radiation would reach space, and chose to absorb 3% of 

the radiation at each of the 15 levels. As can be seen in Fig. 6, this gave 

surprisingly close results. 

To decrease the error associated with the almost flat matrix, the inverse 
of the difference between Dickinson's cooling rates and the initial cooling 
rates was taken. This matrix was also unstable, but it appeared that the 
columns (involving the coefficients scattering to the same level) were 
independent, so we constructed a matrix by averaging by the columns. The 
elements corresponded to ~ 3% absorption at each level and this matrix was 

added to the previous one. But this resulted in a grossly over corrected 

cooling profile and consequently only a fraction of the column- like matrix was 
added to the nearly flat matrix. The resultant matrix yielded cooling rates 
within a few degrees per day of Dickinson's rates, except in the strongly NLTE 
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region. This final Curtis matrix is given in Table 1. This C matrix 
generates a better cooling profile than method 1, but is incapable of 
achieving convergence by method 2. 

An alternate approach would be to use a Curtis matrix calculated 
approximately from first principles and retain K x and K 2 as adjustable 
parameters . 


2. Vertical Constituent Transport in the Mesosphere. 

Another piece of research, partially supported by this grant, was a study 
of vertical constituent transport in the mesosphere in collaboration with NRL 
colleagues (Attachment 1) . Ground-based microwave spectroscopy measurements 
of mesospheric CO and H 2 0 vertical mixing ratio profiles and SME ozone data 
were used to infer vertical mixing rates in the upper mesosphere. The CO and 
H 2 0 data consistently imply vertical eddy diffusion coefficients in the 70-85 
km region of < (1-2) x 10 5 cm 2 s" 1 during spring through summer at 
mid-latitudes. Comparison of SME 0 3 data with model results reinforces the 
conclusion of slow vertical mixing in the upper mesosphere as a consequence of 
the reduced H0 X catalytic loss of odd oxygen. The slow vertical mixing 
deduced in this study is consistent with upper limits obtained from studies of 
the mesospheric heat budget (Apruzese et al., 1984; Strobel et al., 1985) and 

could be construed as evidence for an advectively controlled mesosphere. A 
comparison of the vertical eddy diffusion coefficients for momentum stresses, 
constituent and heat transport suggested that the eddy Prandtl number must of 
order 10. 



5 


References 

Apruzese, J. P. , D. F. Strobel and M. R. Schoeberl, Parameterization of IR 
cooling in a middle atmosphere dynamics model. 2. Non-LTE radiative 
transfer and the globally averaged temperature of the mesosphere and 
lower thermosphere, J. Geophys. Res., 89,, 4917-4926. 

Dickinson, R. E. , Infrared radiative cooling in the mesosphere and lower 
thermosphere, J. Atmos. Terr. Phys., 46, 995-1008, 1984. 

Houghton, J. H. , The Physics of Atmospheres . Cambridge Univ. Press, Cambridge, 
203 pp., 1977. 

Strobel, D. F. , J. P. Apruzese and M. R. Schoeberl, Energy balance constraints 
on gravity wave induced eddy diffusion in the mesosphere and lower 
thermosphere, J. Geophys. Res., 90, 13,067-13,072, 1985. 



6 


Fig. 1 
Fig. 2 

Fig. 3 
Fig. 4 

Fig. 5 

Fig. 6 

Fig. 7 
Fig. 8 

Fig. 9 


Figure Captions 

Temperature T: the temperature data used by Dickinson, in *K. 

Dickinson cooling calculation Q: the desired complex structure, in 
•K/day. 

Approximate source function J: in relative units. 

Modified Curtis Approximation QT: cooling calculated by method 1, 
using the Curtis matrix given in Table 1, in *K/day. 

Error Q-QT: note relatively flat field for all but the upper right 
positions . 

Initial Curtis Approximation QF, using near flat Curtis matrix, in 
•K/day. 

Error Q-QF: worse than Q-QT. 

Iterative Approximation Q5: calculated using near flat matrix and 
only 5 iterations. 

Iterative Approximation QI: calculated using Curtis matrix (Table 1) 
and 5 iterations. 
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